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The physical mechanism of the anthropogenic triggering of large earthquakes on active faults is studied on 
the basis of experimental phenomenology, i.e., that earthquakes occur on active tectonic faults, that crustal 
stress values are those measured in situ and, on active faults, comply to the values of the stress drop measured 
for real earthquakes, that the static friction coefficients are those inferred on faults, and that the effective 
triggering stresses are those inferred for real earthquakes. Deriving the conditions for earthquake nucleation 
as a time-dependent solution of the Tresca-Von Mises criterion applied in the framework of poroelasticity 
yields that active faults can be triggered by fluid overpressures < 0.1 MPa. Comparing this with the 
deviatoric stresses at the depth of crustal hypocenters, which are of the order of 1-10 MPa, we find that 
injecting in the subsoil fluids at the pressures typical of oil and gas production and storage may trigger 
destructive earthquakes on active faults at a few tens of kilometers. Fluid pressure propagates as slow stress 
waves along geometric paths operating in a drained condition and can advance the natural occurrence of 
earthquakes by a substantial amount of time. Furthermore, it is illusory to control earthquake triggering by 
close monitoring of minor "foreshocks", since the induction may occur with a delay up to several years. 

The central role of fluids in seismic faulting is known in all stages of the earthquake cycle, i.e., in the nucleation 
processes, in the dynamic propagation and in the post-seismic evolution. Namely, pore fluid pressure affects 
the whole earthquake process starting from fault nucleation 1,2 , continues through thermally- activated pres- 
surization 37 and mechanical lubrication 8,9 , and finally plays a role in triggering aftershocks 1014 . Interest on the 
first step of the process (fault reactivation and rupture nucleation) has been recently revived by several cases of 
human activities, which are indicted for having induced destructive earthquakes. Davies et al. 15 compiled a list of 
190 possible examples of induced earthquakes, with magnitude spanning between 1.0 and 7.9, connected with 
mining, artificial reservoir impoundment, geothermal operations 16,17 , oil and gas field production and hydraulic 
fracturing (i.e., tracking). 

Most earthquakes induced by anthropogenic activities concern small magnitude events (M < 3.0) located in 
the vicinity of the activities themselves 18,19 . Here we are not interested in such earthquakes, which generally 
constitute more a nuisance than a real danger. Indeed, we are interested in the other type of induced seismicity, 
which regards large events (M > 5.5) on nearby active tectonic faults, at a distance up to a few tens of kilo- 
meters 20,21 . This type of induction is appropriately termed triggered (or activated) seismicity 22 , since human 
activities provide the tiny — but fundamental — input to a system which is independently close to instability. 
A small action produces a large reaction, just as the modest pressure of a finger on a gun trigger releases the large 
explosive energy stored in the propellant. Such a type of induced seismicity has been so far mostly associated with 
the impoundment of artificial water reservoirs, but occasionally also with gas and oil production. 

Triggered seismicity is the most important and dangerous type of induced seismicity, and should not be 
confused with its stablemate, simply called induced seismicity; the latter is commonly associated with drilling 
or hydrofracture and generally implies large local stresses, but small earthquakes. In fact, to induce fresh fracture 
in the bulk rock the applied stresses must be, by definition, equal to those determined by hydrofracture tech- 
niques. However, event size is small because these large stresses are spatially concentrated and can only induce 
fracture on small volumes of rock. On the contrary, triggered seismicity involves large tectonic structures, where 
the stress has been independently accumulating to a near failure conditions by the internal Earth's dynamics 23,24 , 
with the human activity providing only the "last straw". 

Although the unequivocal discrimination between naturally occurring and man-triggered earthquakes is 
difficult 25,26 , there exist some hardly questionable cases. The first documented case of reservoir-triggered seismi- 
city occurred in 1932 in Algeria's Oued Fodda Dam. Another prominent example is provided by the realization of 
the Koyna Dam, India, which was followed in 1967 by the M6.5 Koynanagar earthquake, where 180 people died 



SCIENTIFIC REPORTS | 4 : 6100 | DOI: 1 0. 1 038/srep061 00 



1 



and 1,500 were injured 27 . Other relevant cases are the M6.1 Oroville, 
California, earthquake, also attributed to a reservoir, and the M5.4 
Aswan, Egypt, earthquake in 1981 21 , which occurred 15 years after 
the filling of the Nasser artificial lake. The largest and most recent 
case of reservoir triggered earthquakes is possibly that of the 
Zipingpu Dam in China, which is indicted for triggering the May 
12 2008 M7.9 Wenchuan earthquake, killing some 80,000 people 28-30 . 

Dams are most often at the root of man-triggered seismicity, but 
there are some notable additions. Caused by the injection of pressur- 
ized fluids into a 367-m-deep borehole at the Rocky Mountain 
Arsenal, Colorado was the Denver earthquake sequence in the '60s, 
which culminated in a M5.3 event 31 . The injection of pressurized 
fluids to stimulate the production of an oil field, induced a M5.4 
earthquake in Caviaga, Italy, in 1951 32 , and may have also triggered 
the earthquake sequence of May 20-29 2012 in Emilia, Italy, with two 
events of M6.0 and another five above M5.3, which resulted in the 
death of 27 people as well as a major wound to the gasping Italian 
economy 33,34 . 

The physical problem of quantifying the man-triggering mech- 
anism is a problem still lacking a comprehensive approach, even if 
the basic constitutive equations have been known since quite long 
ago. We work out a solution on the basis of the following phe- 
nomenological ingredients: 

1) Time invariance; since the time scale of earthquake recurrence 
on the same fault segment is 100-1000 years and the evolution 
of the driving mechanism is a geological process 3-4 order of 
magnitudes slower, earthquakes can be taken as repetitions on 
existing weakness zones — the faults. 

2) Stress on faults are close to failure by an unknown extent; the 
stress history on faults depends on the interaction with neigh- 
boring faults and on a variety of factors practically impossible 
to quantify. 

3) Low static friction coefficients; static friction on fault zones is 
much lower than in bulk (undamaged) rocks and in the labor- 
atory (cf. 35 ), and has values ji < 0.3 36-41 . 

4) Very low stress thresholds for triggering; additional fluid pres- 
sures of 0.05 MPa due to artificial water reservoirs have acti- 
vated destructive earthquakes 17 ; similar low values occur in the 
dynamic triggering by nearby earthquakes 42,43 . 

5) Earthquakes triggered with considerable time delay; earth- 
quakes have been triggered with time delays > 10 years 17,21 . 

6) Earthquakes triggered at considerable distance in space; earth- 
quakes have been triggered at up to 30 kilometers 17,21 . 

According to such phenomenological constraints, earthquake 
triggering is likely to occur only on the most favorably oriented 
pre-existent faults, which are at an angle 0 to the largest (in norm) 

principal stress. This angle is given by 0 — - arctg ^— ^ , i.e., according 

to the above-mentioned values of ji, it is comprised in the 37°-45° 
range. Given that the exact [i values are impossible to know before- 
hand, we consider the asymptotic limit fl—>0, which corresponds to 
the maximum shear stress orientation of 45° to the largest principal 
stress. In addition, among the optimally oriented faults, only those in 
close proximity to failure will be triggered. In this framework, we can 
write the failure criterion as 1,44,45 : 

— 2 M — 2 Pflmd ) ( ) 

° eff >°R (2) 

where j.i formally is the static friction coefficient (often written as ji u 
or n s in the literature), a 1 and a 111 are the first and third principal 
stresses, respectively. pnuid is pore fluid pressure, which we will con- 
sider in the range 100-150% of the hydrostatic load (see Numerical 



results and discussion of the failure threshold section). In other 
words, we assume that under normal conditions the fluid is in the 
hydrostatic pressure thanks to a drained percolation-driven circula- 
tion 46 , and we consider the perturbation due to an externally induced 
pore fluid pressure increase. The failure stress a R is a material prop- 
erty, which we assume to be fixed. By keeping /( constant, we impli- 
citly assume that it represents an average value over the whole fault 
surface (i.e. we neglect any possible effect of spatial heterogeneities). 
Moreover, equation (1) also disregards any possible coupling with 
temperature variations and second-order effects 47-49 . 

When fluids are injected into the system, pfl ui d is increased, thus 
increasing a eSi (equation (1)), which may then reach the failure 
condition (equation (2)) 50-52 on a nearby fault which is indepen- 
dently close to failure. We study the problem by first considering a 
time independent model to establish the failure threshold and then 
we consider the time evolution of the process. 

A time independent model 

In reservoirs, water is stored in artificial lakes inducing in the ground 
a fluid pressure equal to the hydrostatic column. Typically, the depth 
h of these artificial basins is some tens of meters and therefore the 
additional overpressure they induce in the subsoil, Pfuad ghisoi the 
order of 0.1 MPa. In a similar manner, pressurized fluids — basically 
water — are injected into ageing gas and oil fields to increase the 
percentage of the extracted oil or gas — the so-called recovery factor - 
or to maintain the productivity over longer periods. In other cases, 
gas (either methane for storage or C0 2 for sequestration and dis- 
posal) is injected in the fluids, resulting in overpressures up to a 
fraction of the ambient fluid pressure. 

While the exact values of pn U jd and its evolution depend on each 
specific case, we consider overpressures up to 50% of the hydrostatic 
pressure: 

Pfluid = (1 + a) phydro =(l + a) p fluid g z (3) 

z being the depth with respect to the ground surface. Such a wide 
range of values of a (a = 0-0.5) can be considered exhaustive of most 
operational configurations. 

Fluid pressures as low as 0.01 MPa have been suggested capable to 
trigger earthquakes 37,53 , while fluid pressures of 0.05 MPa have been 
experimentally found to trigger 17 . We consider a "conservative" 
threshold of 0.1 MPa, a value that can be safely stated as capable of 
earthquake triggering. Comparing this value with the ambient stress 
of ~ 300 MPa at a depth of ~ 10 km, typical of M5.5 crustal earth- 
quake hypocenters, we find that the earthquake triggering stresses 
are comparatively very small (roughly 0.03%), a conclusion consist- 
ent with the values independently inferred for static and dynamic 
aftershock induction 42,43 . 

We take as crustal stress values those determined by in situ mea- 
surements (e.g. 54,55 ; for a comprehensive source see also 56 ), which 
yield that the stresses in the crust are well approximated by a nearly 
lithostatic condition, i.e., by stresses ffy = <5yOV~27 MPa km -1 , and 
deviatoric stresses are of the order of 10 MPa km -1 , or 100 MPa at 
10 km depth. However, the in situ stresses were measured on bulk 
rock and not on faults, even if in some cases not far from them. The 
deviatoric stress values on faults can be directly inferred from the 
stress drop as follows. While the experimental in situ static friction 
coefficients are very low, the dynamic friction coefficients will be 57 
equally low — or lower — implying that earthquakes release virtually 
all deviatoric stresses. Hence, on faults, the latter must be equal to the 
stress drop, which is of the order of 1-10 MPa 58,59 . Using such values, 
the principal stresses on faults remains essentially lithostatic with a 
minor perturbation (cf. equation (4)), i.e. fluid pressure makes faults 
reside always in a state close to failure, in agreement with the crit- 
ically paradigm (for a comprehensive treatment see e.g. 45 ). In other 
words, during an earthquake "cycle" the deviatoric stresses vary only 
within 0.3-3% of the ambient stress, so that any reliable estimate of 
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time-dependent risk of failure would require to know the state of 
stress on the fault with a practically unrealizable accuracy. Hence, we 
must prudentially assume that all faults neighboring the anthro- 
pogenic activity are so close to instability that they will fail when 
a" s increases by 0.1 MPa due to an increment in fluid pressure. 
This condition namely replaces the rupture condition expressed by 
equation (2). 

Numerical results and discussion of the failure 
threshold 

Since our primary goal is to provide a quantitative inference on the 
possible rupture initiation, we merely evaluate if the rupture con- 
dition is met, i.e., if the Tresca-Von Mises criterion is satisfied. We 
use typical in situ stress values 54 " 56 : 



ffi" = l MP km" 



a l + a m = 53 MPa km" 



= 54 MPa km" 



normal faulting 



thrust faulting. 



(4) 



The results are shown in figure 1, where we report the effective stress 
a *as a function of the depth z for a 25% overpressure with respect to 
the hydrostatic pressure. The picture shows 3 that for both compress- 
ive (red curves) and normal faults (blue) earthquake induction is 
possible with static friction coefficients lower than 0.2, with margin- 
ally lower thresholds for the latter. It is also important to note that 
effective stress, and therefore the induction capability, increases lin- 
early with depth, in agreement with the experimental evidence that 
triggered earthquakes tend not to have a superficial source 20,21,34 . 

Figure 2 shows the effective stress as a function of overpressure, 
i.e., of pfl U id normalized to the hydrostatic pressure, at a depth of 5 
kilometers. Earthquake induction appears to be crucially affected by 
both the static friction coefficient and overpressure. The balance is 
subtle, with activation not occurring for faults with a ji larger than 
0.25 at whatever overpressure and occurring only at overpressures 
above 30% for f.i = 0.2. For lower friction coefficients activation is 
possible for all values of overpressure, and failure occurs already at 
0%, i.e., when pore fluid is at the mere hydrostatic pressure, emphas- 
izing the role of precipitation in triggering earthquakes 60 . The activa- 
tion of strike-slip faults is bound to be an intermediate case between 
the normal and compressive cases, and depends on the detail of the 



local conditions. Note that the increase in effective stress reaches 
easily values which are a large fraction of the stress drop — approxi- 
mately 2 MPa in effective stress per 10% increase in fluid pressure 
(see Fig. 2) — implying a possible triggering of all faults, and not only 
of the mature ones, with a substantial advance of the occurrence time. 

Summing up, the time-independent solution candidates fluid 
pressurization as the basic ingredient for earthquake triggering. 
We now proceed to study how the process evolves with time. 

A time-dependent solution 

So far, we have considered a time-independent increase in pore fluid 
pressure, controlled by the overpressure parameter a. The comple- 
mentary time-dependent problem concerns how a time-varying pflum 
at a given spatial source perturbs the fluid pressure field as a function 
of time in the surrounding region. We must solve a coupled elastic- 
hydraulic problem in which fluid migration plays a fundamental role. 
This problem, which can be tackled in the domain of poroleasti- 
city 61,62 , has been already studied by several authors 17,46,63 " 65 in con- 
nection with the seismicity induced by the pore pressure increase due 
to a change in the level of a reservoir lake behind a dam. In a 2-D 
model, the total pore pressure at a depth z, distance x and time r, 
p(r,t), can be expressed as the superposition of the coupled response 
of diffused pore pressure and undrained pore pressure as (ibid.): 



p[r,t)=yp 0 erf 



-po erfc 



2 \fu}i. 



where p 0 is the input pressure at the time t = 0, to the hydraulic 
diffusivity and r = (x 2 + z 2 ) 1 ' 2 the pressure-source to fault distance, 
and y the relative weight of the drained and undrained contributions, 
which in a porous solid under laboratory conditions can be expressed 
as a function of the Skempton coefficient B and of the undrained 

Poisson modulus v„: V = B — , " . 

3 (l-v„) 

However, under real crustal Earth's conditions the effective value 
of y depends on the pattern of faults that rule the draining, which 
occurs essentially through the "easy paths" rather than the bulk 
pores. While knowing beforehand the existence — and the geometry 
— of such a connected pattern is impossible, most authors have 
found that only assuming its existence, i.e. y ~ 0, allows a satisfactory 
fit to the data 10,17,46,63 " 65 . In a similar fashion, while hydraulic diffu- 
sivity of crustal rocks varies experimentally in the laboratory and the 
field over 16 decades 64 , its effective value in the cases of triggered 
seismicity ranges only over 2 decades, from 0.1 to 10 m 2 s _1 46 . 



03 



fc 




z(km) 

Figure 1 | The effective stress <r efi as a function of depth zfor a 25% overpressure with respect to the hydrostatic (i.e., a = 0.25) for different values of the 
static friction coefficient ft . Values of a tfl are computed through equations (1), (3) and (4). Blue curves stand for normal faults and red curves for 
compressive faults. The grey shaded area stands for the region where rupture is not possible. 
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Figure 2 | The effective stress as a function of overpressure for different values of the static friction coefficient // at a depth of 5 kilometers. As in the 

previous figure, blue curves stand for normal faults and red curves for compressive faults. The grey shaded area still denotes the region where rupture is 
not possible. 



We consider here the application to two practical cases, both 
involving an input fluid pressure varying with time, one relative to 
a "short" time scale and another to a "long" time scale. In these cases, 
when the input pore pressure p 0 is time-dependent, the solution can 
be obtained from equation (5) by the principle of superposition at 
different initial times. Figure 3 shows the short time scale case, with 
the input (blue curve) and pore pressures at 1, 2, 4 and 6 km from the 
injection point (red curves), normalized to the input pressure. Its 
main feature — well known in all the cases of triggered earthquakes 
— is the delay with which the pore pressure "wave" propagates. In 
fact, we find that fluid pressure in the crust propagates as a wave: the 



pressure will increase and then die away. Note the difference with a 
simple (and intuitive) steady state solution: in that case pressure will 
increase and then stay at a high value. The peak pressure (marked 
with arrows in Fig. 3) occurs always well after the input pressure 
decrease has begun (in other words, when the pressurization has 
been already terminated); at 1 km it occurs right after the end of 
the full load interval, at 2 km 60 days later, while at 4 km and 6 km it 
occurs nearly 100 and almost 150 days later, when the input pressure 
is below 50% of its maximum. This is important, since it shows how 
early warning procedures to promptly halt the pressurization can 
hardly be effective. Concerning the pore pressures, these reach values 
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Figure 3 | Results pertaining to the fluid diffusion. The input and the resulting pore pressures (blue and red curves, respectively), normalized to the 
former as a function of time at different times and distances r = (xr + z 2 ) 1 ' 2 from the injection point relative to a "short" time scale. Arrows indicate the 
peaks in the resulting pore pressure, which is given by equation (5). The input pressure increases from 0 to its maximum according to a linear ramp 42 
days long, it stays constant for 26 days and goes back to 0 according to a linear ramp lasting 1 84 days. The horizontal portion of this input signal (i.e., where 
it reaches its maximum value) is associated with the values given by equation (3), and thus it is ultimately related to a specific value of the overpressure 
parameter i. More complicated input time functions can be handled, but the present choice, which describes realistic operational conditions, makes it 
possible to capture the most important features of the problem. A value of y = 0.01 was assumed together with a hydraulic diffusivity oj = 10 m 2 s _1 . 
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Figure 4 | The same as Fig. 3, but relative to a "long" time scale. It concerns an input pressure increasing from 0 to its maximum according to a linear 
ramp 2.8 years long, followed by a constant load for 1.7 years and going back to 0 according to a linear ramp lasting 12 years. A value of y = 0.0001 was 
assumed together with a hydraulic diffusivity m = 0.1 m 2 s -1 . 



larger than 20% of the maximum input pressure at distances less than 
2 km, but are respectively one and two orders of magnitude smaller 
at 4 and 6 km. The pore pressure curve at 6 km from the injection 
point shows initially the predominance of the undrained elastic term, 
with an increase almost simultaneous with the load, and then, after 
more than 100 days, the delayed arrival and predominance of the 
hydraulic pressure wave. 

Figure 4 shows the "long" time scale case, with the input and pore 
pressures computed at 3, 6, 12 and 18 km from the injection point, 
still normalized by the input pressure. Again, the picture is domi- 
nated by the delayed propagation of the pore pressure wave, which 
reaches its peak value 8 and 10 years after the input at 12 km and 
18 km, respectively, attaining pressure values which are within some 
percent of the input. Remarkably, the latter values are only appar- 
ently small, if one recalls that hydrostatic pore fluid pressure at 
hypocentral depths is 3 orders of magnitude larger than the stresses 
found to trigger earthquakes near reservoirs (see A time independent 
model section). 

Discussion and conclusions 

Considering the case of induced seismicity relative to large destruct- 
ive earthquakes, we found that they are due to small fluid overpres- 
sures which trigger nearby tectonic faults. Pulling the trigger of 
loaded gun is an action minute per se, but decisive for shooting. 
The same earthquakes would occur naturally, but at a (possibly 
much) later time, due to an increase in the shear stress induced by 
tectonic loading, to stress perturbations from slip on neighboring 
faults, to the fluid overpressure of a nearby earthquake, or the 
dynamic load imposed by the waves of the same. The time required 
to naturally achieve the failure depends on a variety of factors, start- 
ing from the load rate and the residual stress left over by previous 
earthquakes. Since the "recurrence" time of earthquakes on the same 
fault segment is of the order of 10 3 years 66 68 and since the stress drop 
is of the order of 1-10 MPa 58,59 ' 69 ' 70 , the threshold we used places the 
lower limit of anthropogenic advance in the time of occurrence at 
some decades, but the calculated values of the effective stresses yield 
that the clock advance (see, e.g., Fig. 2 oP 4 and references cited 
therein) can be very substantial. 



The effect of this anthropogenic perturbation is like anticipating 
the death of a living organism to a sooner date than that which would 
naturally occur; as such, it carries some potential societal respons- 
ibilities. Another important issue emerging from the above results is 
that there is little hope to control triggered seismicity through an 
early warning based on intensive monitoring: the pore pressure wave 
propagates with such a delay that no matter how fast one reacts to the 
onset of minor precursory seismicity, he will be too late (see A time- 
dependent solution section). In fact, the fluid pressure pulse propa- 
gates independently of what is done at the source. Actually, it may 
still be propagating even after the whole apparatus that produced it 
has been not only shut off, but also dismantled. 

To recap, as pointed out by 71 and reiterated by 31,72 75 and many 
others, the fact that fluid injection into the ground may trigger tec- 
tonic earthquakes has been known for quite long 76,77 , although the 
combined lack of a thorough comprehension of the basic mechan- 
isms has possibly prevented considering them always with the due 
caution. Indeed, a general policy which maximizes the oil and gas 
recovery and storage and concurrently minimizes the risks of fluid 
injection has still to be defined. Any project should primarily include 
a detailed mapping of nearby active faults, carefully evaluating the 
seismic risk and — if the case — proceeding to secure the region prior 
to any operation. 
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